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Abstract 

A first-ever 2-dimensional celestial map of primary cosmic-ray flux was ob- 
tained from 2.10 x 10 8 cosmic-ray muons accumulated in 1662.0 days of Super- 
Kamiokande. The celestial map indicates an (0.104 ± 0.020)% excess region 
in the constellation of Taurus and a —(0.094 ± 0.014)% deficit region toward 
Virgo. Interpretations of this anisotropy are discussed. 



* Talk at "Les Rencontres de Physique de la Vallee d'Aosta (La Thuilc 
2006)", La Thuile, Aosta Valley, Italy, March 5-11, 2006. 

t The talk is based on G.Guillian et al. (Super-Kamiokande collaboration), 
submitted to Phys.Rev.D, |astro -ph/0508468. 

* The PowerPoint file used in the talk can be downloaded from 
http : / / www-nu . kek . j p/ "oyama/LaThuile . oyama . ppt 



1 Super-Kamiokande detector and cosmic-ray muon data 

Super-Kamiokande (SK) is a large imaging water Cherenkov detector located 
at ^2400 m.w.e. underground in the Kamioka mine, Japan. The geographical 
coordinates are 36.43°N latitude and 137.31°E longitude. Fifty ktons of water 
in a cylindrical tank is viewed by 11146 20-inch</> photomultipliers. 

The main purpose of the SK experiment is neutrino physics. In fact, SK 
has reported many successful results on atmospheric neutrinos and on solar 
neutrinos. For recent results on neutrino physics as well as the present status 
of the SK detector, see Koshio.D 

The SK detector records cosmic-ray muons with an average rate of ~1.77 Hz. 
Because of more than a 2400 m.w.e. rock overburden, muons with energy larger 
than ~1 TeV at the ground level can reach the SK detector. The median en- 
ergy of parent cosmic-ray primary protons (and heavier nuclei) for 1 TeV muon 
is ~10 TeV. 

Cosmic-ray muons between June 1, 1996 and May 31, 2001 were used in 
the following reported analysis. The detector live time was 1662.0 days, which 
corresponds to a 91.0% live time fraction. The number of cosmic-ray muons 
during this period was 2.54 x 10 s from 1000 m 2 ^1200 m 2 of detection area. 

Muon track reconstructions were performed with the standard muon fit 
algorithm, which was developed to examine the spatial correlation with spal- 
lation products in solar neutrino analysis. ^ In order to maintain an angular 
resolution within 2°, muons were required to have track length in the detector 
greater than 10 m and be downward-going. The total number of muon events 
after these cuts was 2.10 x 10 s , corresponding to an efficiency of 82.6%. 

2 Data analysis and results 

The muon event rate in the horizontal coordinate is shown in Fig^ The rate 
is almost constant and the time variation is less than 1%. This distribution 
merely reflects the shape of the mountain above the SK detector. For example, 
the muon flux from the south is larger because the rock overburden is small in 
the south direction. 

With the rotation of the Earth, a fixed direction in the horizontal coordi- 
nate moves on the celestial sphere. Therefore, the time variation of muon flux 
can be interpreted as the anisotropy of primary cosmic-ray flux in the celes- 
tial coordinate. ^ A fixed direction in the horizontal coordinate travels on a 
constant declination, and returns to the same right ascension after one sidereal 
day. The muon flux from a given celestial position can be directly compared 
with the average flux for the same declination. 

Since 360° of right ascension is viewed in one sidereal day, the right- 
ascension distribution is equivalent to the time variation of one sidereal day 



period. The cosmic-ray muon flux may have other time variations irrelevant of 
the celestial anisotropy, for example, a change of the upper atmospheric tem- 
perature,® or the orbital motion of the Earth around the Sun. An interference 
of one day variation and one-year variation may produce a fake one sidereal 
day variation. Those background time variations are carefully examined and 
removed to extract ~ 0.1% of the real primary cosmic-ray anisotropy. For more 
details, see G.Guillian et al. ® 

The deviations of the muon flux from the average for the same declination 
are shown in Fig. [2] The units are amplitude in Fig. HJa) and significance in 
Fig. Gib). Obviously, an excess is found around a p=> 90° and an deficit around 
a re 200°. (The excess and deficit around <5>70° and S< - 40° in Fig^a) are 
due to poor statistics, as can be recognized from Fig.^b).) 

To evaluate the excess and deficit more quantitatively, conical angular 
windows are defined with the central position in the celestial coordinates (a, S) 
and the angular radius, A8. If the number of muon events in the angular- 
window is larger or smaller than the average by 4 standard deviations (which 




Figure 1: Cosmic-ray muon rate in the horizontal coordinate. The units are 
day _1 m _2 sr _1 . The dot curves indicate contours of constant declination, while 
the arrows indicate the apparent motion of stars with the rotation of the Earth. 
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Figure 2: Primary cosmic-ray flux in the celestial coordinate. Deviations 
from the average value for the same declinations are shown. The units are 
(a)amplitude (from —0.5% to 0.5%) and (b)significance (from — 3er to 3a). 
The Taurus excess is shown by the red solid line and the Virgo deficit is shown 
by the blue solid line. 



Table 1: Amplitude, center of the conical angular windows in the celestial 
coordinate, angular radius of the window, the chance probability of finding the 
excess or deficit are listed. Small chance probabilities might occur somewhere 
on the map because all positions on the celestial sphere are surveyed. Chance 
probabilities considering such a "trial factor" are also listed in the last row. 



Name 


Taurus Excess 


Virgo deficit 


Amplitude 


(1.04 ±0.20) x 10- 3 


-(0.94 ±0.14) x 10" a 


Center (a, S) 


(75° ±7°, -5° ±9°) 


(205° ±7°, 5° ±10°) 


Angular radius (A9) 


39° ± 7° 


54° ± 7° 


Chance probability 


2.0 x 10~ 7 


2.1 x 10" 11 


(trial factor is considered) 


5.1 x 10~ 6 


7.0 x 10~ n 



corresponds to chance probability of 6.3 x 10 -5 ), the angular window is defined 
as the excess window or the deficit window. The celestial position (a, S) and 
the angular radius (A9) are adjusted to maximize the statistical significance. 

By this method, one significant excess and one significant deficit are 
found. From the constellation of their directions, they are named the Tau- 
rus excess and the Virgo deficit. Summary of the Taurus excess and the Virgo 
deficit are listed in Table 1. The positions of the Taurus excess and the Virgo 
deficit are also shown in Fig. [3 

3 Comparison with other experiments 

Fig. |21 is the first celestial map of cosmic-ray primaries obtained from under- 
ground muon data. However, there are three similar celestial maps from other 
experiments, even though they are not published in any refereed papers. Two 
of them are from 7-ray observatories: Tibet air shower 7 observatory E» and 
Milagro TeV-7 observatory. 13 Note that their primary particles include not 
only protons, but also 7-rays, because of their poor proton/7 separation capa- 
bility. The other is a celestial map from the 1MB proton decay experiment.^ 
In the 1MB map, only excess regions are plotted. Results from the 3 experi- 
ments are shown in Fig. [21 The trends of 3 celestial maps well agree with the 
results from SK. 

In addition to three celestial maps, there were many one-dimensional 
results from underground cosmic-ray muon observatories. Most of the exper- 
iments use very simple detectors, such as 2 or 3 layers of plastic scintillators. 
They count cosmic-ray muon rate with coincidence of the plastic scintillator 
layers. All cosmic-ray muons are assumed to arrive from the zenith, and right- 
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Figure 3: Primary cosmic-ray flux distribution from 3 experiments. They 
are Tibet air-shower 7 observatory (top), Milagro TeV-7 observatory (mid- 
dle), and 1MB proton decay experiment (bottom). © The center of the 
Taurus excess by SK is indicated by a star, and the center of Virgo deficit is 
indicated by a triangle. 



ascension distributions are fitted with first harmonics. The declination distri- 
bution cannot be analyzed. 

Exactly the same analysis method was applied to the SK data to examine 
the consistency. The right-ascension distribution is shown in Fig. 0] The 
amplitude and the phase of the first harmonics were obtained to be (5.3 ± 
1.2) x 10~ 4 and 40° ± 14°. Results of the analysis are plotted together with 
other underground muon experiments and some air shower array experiments 
in Fig. [SJ The agreement with other experiments is excellent. Especially, the 
phases of most experiments range between 0° and 90°. 




Figure 4: Cosmic-ray muon rate as a function of the right ascension in Super- 
Kamiokande. The average muon rate is normalized to be 1. It is assumed 
that all muons come from the zenith. The solid curve is the best fit of the 
first two harmonic functions. The dashed curve is the first two harmonics 
after subtracting the atmospheric contribution (See Guillian et al. -'). The 
amplitude and phase of the first harmonics are (5.3 ±1.2) x 10 -4 and 40° ± 14°. 
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Figure 5: First-harmonic fit of right ascension distributions by various cosmic 
ray experiments. The amplitude (top) and phase (bottom) are plotted as 
a function of the primary energies. The circles are for underground muon 
experiments and squares are for extensive air shower arrays. The filled circle is 
for Super-Kamiokande. Data references are as follows: Bo:Bolivia(vertical),l3 
Mi:Misato(vertical), CHJ Bu:Budapest, Hob:Hobart(vertical), 
Ya:Yakutsk, E3) LoV:London(vertical), EH So:Socomo(vertical), © 
Sa:Sakashita(vertical) ,EJ LoS:London(south) , E3> Li:Liapootah(vertical) , 
Ma:Matsushiro(vertical), Ell Ot:Ottawa(south), E3> P :Poatina(vertical), 
Ho:Hong Kong, Ut:Utah, ^ BaS:Baksan(south), 

Kam:Kamiokande, E3> Mac:MACRO, 12} Tib:Tibet(vertical), 
Ba:Baksan air shower, No:Mt. Norikura, Ea:EAS-TOP, ^ 

Pe:Peak Musala. ESJ "(vertical)" means that the upper plastic scintilla- 
tor layers are placed exactly above the bottom layers and the coincidence is 
sensitive to muons from the zenith, "(south)" means that the upper layers are 
placed rather south of the bottom layers and muons from south direction are 
selectively counted. 



4 Can protons be used in astronomy? 



Before interpretations of the SK cosmic-ray anisotropy, trajectories of protons 
in the galactic magnetic field must be addressed. The travel directions of 
protons are bent by the galactic magnetic field in the Milky Way, which is 
known to be ~ 3 x 10 _ 10 Tesla. If the direction of the magnetic field is vertical to 
the proton direction, the radius of curvature for 10 TeV protons is ~ 3x 10 -3 pc. 

Since the radius of the solar system is ~ 2 x 1CP 4 pc, 10 TeV protons 
keep their directions from outside of the solar system. On the other hand, 
since the radius of the Milky Way galaxy is ~ 20000 pc, protons may loose 
their directions on the scale of the galaxy. 

However, if the magnetic field is not vertical to the proton direction, 
the trajectories of protons in a uniform magnetic field become spiral. The 
momentum component parallel to the magnetic field remains after a long travel 
distance. Since the galactic magnetic field is thought to be uniform on the order 
of ,>300 pc, protons may keep their directions within this scale. The actual 
reach of the directional astronomy by protons is unknown. 

5 Excess/deficit and Milky Way galaxy 

Directional correlations of Taurus excess and Virgo deficit with the Milky Way 
galaxy is of great interest. Schematic illustrations of Milky Way galaxy are 
shown in Fig|J3 Milky Way is a spiral galaxy with 20000 pc radius and ,> 200 pc 
thickness. The solar system is located about 10000 pc away from the center of 
the galaxy. It is in the inside of the Orion arm and about 20 pc away from the 
galactic plane, as shown in FigEfbottom) . 

The Taurus excess is toward the center of the Orion arm, and the Virgo 
deficit is toward the opposite to the galactic plane. Accordingly, primary cosmic 
ray flux have a positive correlation with density of nearby stars around the 
Orion arm. 



6 Compton-Getting effect 

Assume that "cosmic-ray rest system" exists, in which the cosmic-ray flux is 
isotropic. If an observer is moving in this rest system, the cosmic-ray flux 
from the forward direction becomes larger. The flux distribution ($(0)) shows 
a dipole structure, which is written as $(0) oc 1 + A cos (9, where 9 is the 
angle between the direction of the observer's motion and the direction of the 
cosmic-ray flux. Such an anisotropy is called the Compton-Getting effect. 
The velocity of the observer (v) is proportional to A. If v is 100 km/s, A is 
1.6 x 10-3. 







Central bulge 
/ , Galactic nucleus 






Sun 
1 1 








1— 1 

Disk 


• 

^> Globular dusters 

• 




Hale 


" ■ : • 





Orion arm 



Virgo deficit 



^lOOpcl^ggf J ~20pc 



Taurus excess 



~1000pc 



Figure 6: Top view (top) and side view (middle) of the Milky Way galaxy. 
The position of the solar system is shown by red circles. A cross-sectional view 
of the Orion arm around the Earth is also shown (bottom). The Orion arm 
is ^1000 pc in width and >200 pc in thickness. The solar system is inside of 
the Orion arm and ~20 pc away from the center of the Galactic plane. The 
direction of the Taurus excess and the Virgo deficit are also shown. 



If the Taurus excess (1.04xlCn 3 ) and the Virgo deficit (— 0.94xl0~ 3 ) were 
in opposite directions, it might be explained by the Compton-Getting effect 
of v — 50 ~ lOOkm/s. However, the angular difference between the Taurus 
excess and the Virgo deficit is about 130°. The Taurus- Virgo pair is difficult 
to be explained by the Compton-Getting effect. Accordingly, a clear Compton- 
Getting effect is absent in the SK celestial map. Although it is difficult to set 
an upper limit on the relative velocity because there exist excess and deficit 
irrelevant to Compton-Getting effect, it would be safe enough to conclude that 
the relative velocity is less than several ten km/s. 

The relative velocity between the solar system and the Galactic center 
is about 200 km/s. The velocity between the solar system and the microwave 

background is about 400 km/s. The velocity between the Milky Way and 

the Great Attractor is about 600 km/s.ES' The upper limit, several ten km/s, 
is much smaller than those numbers. The cosmic-ray rest system is not to- 
gether with the Galactic Center nor the microwave background nor the Great 
Attractor, but together with our motion. 

Because of the principal of the SK data analysis, two possibilities cannot 
be excluded: the Compton-Getting effect is canceled with some other excess 
or deficit, and the direction of the observer's motion is toward S ~ 90° or 
S ~ -90°. 

7 Crab pulsar 

One strong interest concerns the correlation of the Taurus excess with the Crab 
pulsar. This provides a clue to examine whether high-energy cosmic rays are 
accelerated by supernovae or not, which is the most fundamental problem in 
cosmic-ray physics. 

The Crab pulser^J j s a neutron star in the Crab Nebula, which is one 
of the closest and newly exploded supernova remnants. The distance from the 
solar system is about 2000 pc, and the supernova explosion was in 1054. The 
celestial position is (a, S) — (83.63°, 22.02°). It is in the constellation Taurus 
and also within the Taurus excess, but is deviated from the center of the Taurus 
excess by 28°. 

The total energy release from the Crab pulsar is calculated from the 
spin-down of the pulsar, and is 4.5 x 10 38 erg-s _1 . If it is assumed that all 
energy release goes to the acceleration of protons up to 10 TeV and the emis- 
sion of protons is isotropic, the proton flux at the Earth is calculated to be 
- 0.6 x 10 -7 cm- 2 s -1 

On the other hand, the Taurus excess observed in SK is converted to the 
primary proton flux at the surface of the Earth using the observation period 
and the detection area. The flux is obtained to be ~ 1.8 x 10 _7 cm~ 2 s _1 



From a comparison of these two numbers, the Taurus excess cannot be 
explained by the proton flux accelerated by Crab pulsar by a factor of ~ 3. Note 
that two extremely optimistic assumption were implicitly made in calculating 
the expected flux from Crab pulsar; all energy releases are provided to the 
acceleration of protons up to 10 TeV, and protons travel straight to the Earth. 

8 Summary 

The first-ever celestial map of primary cosmic rays (> 10 TeV) was obtained 
from 2.10x 10 8 cosmic-ray muons accumulated in 1662.0 days of Super-Kamiokande 
between June 1, 1996 and May 31, 2001. In the celestial map, one excess and 
one deficit are found. They are (1.04 ±0.20) x 10 -3 excess from Taurus (Taurus 
excess) and —(0.94 ±0.14) x 10~ 3 deficit from Virgo (Virgo deficit). Both of 
them are statistically significant. Their directions agree well with the density 
of nearby stars around the Orion arm. A clear Compton-Getting effect is not 
found, and the cosmic-ray rest system is together with our motion. The Taurus 
excess is difficult to be explained by the Crab pulsar. 

In 1987, Kamiokande started new astronomy beyond "lights". In 2005, 
Super-Kamiokande started new astronomy beyond "neutral particles" . 

-Note Added- 

After the submission of the first draft, » it was pointed out that an in- 
terpretation as Compton-Getting effect is not impossible. Although the angle 
between the centers of the Taurus excess and the Virgo deficit is 130°, agree- 
ment with the dipole structure is fair (not good) because of the quite large 
angular radius of the window (39° for the Taurus excess and 54° for the Virgo 
deficit). Even if the Taurus- Virgo pair is due to the Compton-Getting effect, 
the relative velocity is about 50 km/s. This does not change the discussion 
about the comparison with other relative velocities. 
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